两种高分辨率数值方法计算波前有非均匀流的变截面管内激波传播规律的比较

本文采用两种对间断解具有高分辨率的数值方法基于推广Riemann问题解的二阶Godunov型有限差分法和分裂算子的随机选取法,计算了微波在波前有非均匀定常流的一维变截面管道中的传播和波后流场特征,得到一致结果,用数值模拟方法揭示了这类运动的一些特殊规律。对比两种方法的计算过程和结果,可以看出,二阶Godunov型方法明显优于随机选取法。     

Numerical Simulation of Plasma Dynamics in a Nonuniform Magnetic Field

An efficient algorithm is proposed enabling numerical simulations of plasma dynamics in a nonuniform magnetic field. The present numerical data are in good agreement with experimental data obtained in a GOL-3 setup and with previous simulations. The experimentally observed effect of fast transfer of energy to ions is confirmed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 1, pp. 35–45, January–February, 2006.     

基于CUDA的有限元矩阵并行装配算法研究

构建航天飞行器的结构有限元模型是准确模拟飞行仿真、完成飞行器在轨飞行阶段结构故障监测和诊断的基础。采用细长体飞行器简化梁模型，提出新的基于CUDA（Compute Unified Device Architecture）的有限元单元刚度矩阵生成和总刚度矩阵组装算法。依据梁单元矩阵的对称性，结合GPU硬件架构提出并行生成算法并进行改进。为有效减少装配时间，在装配过程中采用着色算法，提出了基于GPU（Graphics Processing Unit）共享内存的非零项组装策略，通过在不同计算平台下算例对比，验证了新算法的快速性。数值算例表明，本文算法的求解效率较高，针对一定计算规模内的模型可满足快速计算与诊断的实时性要求。     

Nonlinear stability analysis of miscible displacements

A simulator for three-dimensional horizontal miscible displacements in porous media is developed. Using this simulator, we examine the initiation and development of instabilities, viscous fingers and gravity tongues.With the only perturbations to the system being truncation and round-off errors, a density ratio (the ratio of the density of the displacing fluid to that of the displaced fluid) different from one is responsible for the initiation of the instabilities, and an unfavorable mobility ratio (the ratio of the viscosity of the displaced fluid to that of the displacing fluid) is responsible for the growth of the instabilities.     

三维有限体中平片裂纹的超奇异积分方程方法—Ⅱ.数值方法

超奇异积分方程方法的理论分析已在本文的第Ｉ部分中给出，这一部分是经的数值方法，及用此方法求解的若干典型的平片裂纹问题。     

氢氧混合气体爆轰波的真实化学反应模型数值模拟

采用高精度的ENO格式和基于基元化学反应的真实化学反应模型求解氢氧混合气体一维爆轰波的精细结构。采用直接起爆方法得到稳定传播的爆轰波 ,计算的爆轰波阵面参数和实验相当符合。对爆轰波反应区化学反应的研究表明 ,参与反应的不同组分具有不同类型的变化特征。网格尺寸影响的研究表明 ,计算结果的精度随着网格尺寸的增加而增加 ,并能保持较好的收敛性。移动网格研究结果表明 ,网格运动速度和爆轰速度接近时 ,两者的相互作用对计算结果产生一定影响。     

土中爆炸地冲击能量分布研究

利用数值模拟方法研究了土中装药不同埋设深度爆炸能量分布问题,通过试验得到了耦合系数数据。给出了自由场条件下,封闭爆炸最小比例埋深为2.0 m/kg1/3,这与美军设计规范TM5-855-1中的数据0.56 m/kg1/3相差较大,分析了产生差异的可能原因。并指出本文中给出的冲量型耦合系数与TM5-855-1中的应力型耦合系数是有所区别的,应该注意其使用范围和对象。     

NUMERICAL SIMULATION OF OIL MIGRATION-ACCUMULATION OF MULTILAYER AND ITS APPLICATION

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｏｉｌｆｏｒｍａｔｉｏｎｉｎｓｅｄｉｍｅｎｔｂａｓｉｎｓ ,ｉｔｓｄｉｓｐｌａｃｅｍｅｎｔ,ｔｒａｎｓｐｏｒｔａｎｄａｃｃｕｍｕｌａｔｉｏｎ ,ａｎｄｔｈｅｆｉｎａｌｆｏｒｍａｔｉｏｎｏｆｏｉｌｄｅｐｏｓｉｔｓｈａｖｅｂｅｅｎｏｎｅｏｆｔｈｅｋｅｙｐｒｏｂｌｅｍｓｉｎｔｈｅｅｘｐｌｏｒａｔｉｏｎｏｆｏｉｌ＿ｇａｓｒｅｓｏｕｒｃｅｓ.Ｈｏｗｈａｓｏｉｌｂｅｅｎａｃｃｕｍｕｌａｔｅｄｉｎｔｈｅｐｒｅｓｅｎｔｌｏｏｐａｃｃｏｒｄｉｎｇｔｏｔｈｅｍｅｃｈａｎｉｃｓｏｆｉｍｍｉｓｃ…     

爆炸作用数值模拟

采用简单的炸药化学反应率方程，利用二维流体弹塑性流动程序，对一个变直径铜壳中装药的爆炸作用做了数值模拟。     

Finite element analysis for the unsteady nearshore circulation due to wave-current interaction (I)—numerical model

In this paper, a numerical model for predicting the unsteady nearshore circulation due to wave-current interaction was proposed. In addition to the traditional continuity, momentum and energy equations, the dispersion and refraction relations were included in the governing equations. Moreover, the effects of lateral shears, wind, radiation and bottom stresses were analysed in the governing equations. Therefore, we expect that this model may more completely and exactly reflect the law of wave-current interaction. In part (II) we will adopt the selective lumping two-step explicit finite element method to solve the model, and some examples will be presented.